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Abstract 

We describe microscopic theory for the quantum transport through 
finite interacting systems connected to noninteracting leads. It can 
be apphed to smah systems such as quantum dots, quantum wires, 
atomic chain, molecule, and so forth. The Keldysh formalism is in- 
troduced to study the nonlinear current-voltage characteristics, and 
general properties of the nonequilibrium Green's functions are pro- 
vided. We apply the formulated to an out-of-equilibrium Anderson 
model that has been used widely for the Kondo effect in quantum 
dots. In the linear-response regime the Kubo formalism still has an 
advantage, because it relates the transport coefficients directly to the 
correlation functions defined with respect to the thermal equilibrium, 
and it has no ambiguities about a profile of the inner electric field. 
We introduce a many-body transmission coefficient, by which the dc 
conductance can be expressed in a Landauer-type form, quite gener- 
ally. We also discuss transport properties of the Tomonaga-Luttinger 
liquid, which has also been an active field of research in this decade. 
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1 Introduction 



Quantum transport through finite interacting-electron systems has been stud- 
ied extensively in this decade. For instance, the Coulomb blockade and var- 
ious effects named after Kondo, Aharanov-Bohm, Fano, Josephson, etc. in 
quantum dot systems have been a very active field of research. Furthermore, 
the realization of non-Fermi liquid systems such as the Tomonaga-Luttinger 
liquid in quantum wires and multi-channel Kondo behavior in some novel 
systems have also been investigated by a number theorists. 

To study transport properties of correlated electron systems, theoretical 
approaches that can treat correctly both the interaction and quantum inter- 
ference effects are required. The Keldysh Green's function approach is one 
of such methods [Ill2ll3lllElinilZllHlini- Specifically, the formulation for the 
nonlinear current-voltage profile by Caroli et al has been applied widely to 
the quantum transport phenomena. In this report, we describe the outline 
of the Keldysh formalism in Sec. |2l and then apply it to a single Anderson 
impurity, which is a standard model of quantum dots in the Kondo regime 
in Sec. El 

When a finite sample is connected to reservoirs that can be approximated 
by free-electron systems with continuous energy spectrums, the low-energy 
eigenstates of whole the system including the attached reservoirs are deter- 
mined coherently. Thus, to understand the low-temperature properties, the 
information about the low-lying energy states of the whole system is required. 
The local Fermi-liquid theory [TIH ITT] , which was originally introduced for 
the Kondo systems |12j, is also applicable to the transport properties in wide 
classes of the interacting-electron systems at low temperatures. In Sec. HI 
we reformulate the transport theory for the interacting systems connected 
to noninteracting leads based on the Kubo formalism. The dc conductance 
can be written in a Landauer-type form with a many-body transmission co- 
efficient determined by a three-point correlation function. We also provide 
a brief introduction to Tomonaga-Luttinger model in Sec. El to take a quick 
look at the transport properties of a typical interacting-electron system in 
one dimension. 
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Left lead Sample Right lead 




Figure 1: Schematic picture of the system. 



2 Keldysh formalism for quantum transport 
2.1 Thermal equilibrium 

We start with a system that consists of three regions: a finite central region 
(C) and two reservoirs on the left (L) and the right {R). The central region 
consists of resonant levels, and the interaction Uj^j-^-j^ji is switched on 
only for the electrons in this region. We assume that each of the reservoirs 
is infinitely large and has a continuous energy spectrum. The central region 
and reservoirs are connected with the mixing matrix elements Vj^ and f^, as 
illustrated in Fig. ^ The complete Hamiltonian is given by 

'^tot = T^L + T^R + T^c + + T^mix ■ (1) 

^L = -y: t% cic^, , hr = -Y^^z 4s-. ' (2) 

a (T 

ij&c ^ {i}ec 

^mix = - I] ^^L [c^C-la + H.C. \-Y^R [cjv+la C^Va + H.C. ] . (4) 
cr (7 

Here, c]^ (c^-^) creates (destroys) an electron with spin a at site j, and [i 
is the chemical potential. Also, t^, t^, and t,^ are the intra-region hopping 
matrix elements in each of the three regions L, i?, and C, respectively. The 
labels 1, 2, . . ., are assigned to the sites in the central region. Specifically, 
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Figure 2: Examples for the profile of the electrostatic potential ^ci'i) for 
(a) an insulating sample, and (b) a metal sample, where eV = — ^r. 



the label 1 (A^) is assigned to the site at the interface on the left (right), and 
the label (A^ + 1) is assigned to the site at the reservoir-side of the left 
(right) interface. We will be using units h = 1 unless otherwise noted. 
The density matrix for the equilibrium state peq is given by 

ctCia , Nc = 4 , Nn = 4aCia ■ (6) 

ieL,(T jeC,(T i£R,a 

Therefore, the Hamiltonian Tilat aiid a single chemical potential /j, determine 
the statistical weight in thermal equilibrium. 

2.2 Statistical weight for nonequilibrium steady states 

When the voltage V is applied, the contribution of the electrostatic potential 
has to be included into T^tot^ as 

'Htot = "^tot + "^'ext , (7) 
Vext = ^lNl + ^rNr+ J2 "^Ci^ ■ (8) 

Here $l and $r are the potentials for the lead at L and R, respectively, 
and the applied bias voltage corresponds to eV = — ^r- To determine 
the potential profile in the central region ^^(z), the energy of the electric 
field should also be included into the Hamiltonian eq. (j21), and it should be 
determined self-consistently. However for simplicity, we assume that $c(«) 
is a given function. In Fig. |2l two typical profiles are illustrated. For an 
insulating sample, there must be a finite electric field in the central region 



Pcq 
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and the potential shows approximately a linear z-dependence as that in the 
panel (a). In an opposite case, for a metaUic sample, the electric field vanishes 
in the central region, and the potential profile will become the one as shown 
in the panel (b). Realistic situations seem to be in between these two extreme 
cases. 

In contrast to the thermal equilibrium, one cannot write down the density 
matrix that describes the nonequilibrium statistical weight simply by using 
a single chemical potential, 

p Q-P{T<tot-KNL+Nc+NR)} I rj^ ^-p{H,,^-t_i{NL+Nc+NR)} ^ ^g^ 

because this statistical weight describes the situation after the electrons have 
already been redistributed to gain the electrostatic potential energy. One 
possible statistical weight that describes a nonequilibrium steady state was 
introduced by Caroli et al jl]. 

In the formulation of Caroli et al, the coupling to the leads Timix and 
the interaction Ti^ in the sample region are switched on adiabatically by 
separating the total Hamiltonian in the form 



n 



it) 

n2{t) 



Hi + n2{t) 



R 



Hi: 



R 



n mix 



(10) 

(11) 

(12) 
(13) 

(14) 



Here, 5 = 0^ is an positive infinitesimal. Because Tii has a quadratic form, 
it is possible to use the Wick theorem in the perturbation expansion with 
respect to Ti.2- At t = — oo the two reservoirs and the impurity are isolated, 
so that the different chemical potentials fiL, fiR, and fic can be introduced 
into the three regions in the initial condition. The time evolution of the 
density matrix is determined by the equation 



^tot(t) , P{t) 



(15) 



-iHit 



as 



The formal solution of this equation can be obtained, by using the interaction 
representation p{t) = e''^^^ p{t) e''^^^ and 7^2(t) = e^^^*H2(t) e 

p{t) = U{t,to)pito)U{to,t) , (16) 
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U (t, to) = T exp 
U{tQ,t) = Texp 



to 



dt' rtiit') 



to 



dt' H2{t') 



(17) 
(18) 



Here, T denotes the operator for the chronological time order, and T is the 
anti-time-ordering operator. Note that eq. (|18p is the Hermite conjugatae of 
eq. p7|) . Caroli et a/ have assumed that the initial condition at to ~oo is 
given by 



p(-cx)) 



g-/3{Wi;L-MLAfL} g-/3{Wi;C-McAfc} P{Hl;R- r} 



Tr 



'P{'Hl■,L-^^LNL} p-PiHi.c-t^cNc} e-l^{'^uR-^^RNR} 



(19) 



Namely, at to ~^ —oo, each system is in a thermal equilibrium state with the 
chemical potential iil,r,c- Here fJ^L — fJ^R = — = . 



2.3 Perturbation expansion along the Keldysh contour 

The perturbed part 'H2{t) is switched fully on at t = 0. Therefore the ex- 
pectation value of physical quantities are defined with respect to the density 
matrix at t = 0, 

{O) ^ Tr[p(0)O5] 

= Tr[p(0)O5] = Tr[p(-oo)f/(-oo,0)O5?7(0,-oo)] , (20) 

where Os is a Schrodinger operator, and [p(0), Htot(O)] = for the sta- 
tionary states. Equation can be rewritten by using a property of the 
time-evolution operator, [/(— cx),0) = f/(— oo, +oo) t/(-|-oo, 0), as 

(C) = ( f/(-oo, +oo) f/(+oo, 0) Os f/(0, -oo) )o , (21) 

where (■ ■ ■)q = Tr [p(— oo) • ■ ■]. The stream of the time in eq. ()21|) can be 
mapped onto a loop shown in Fig. IHl starting from t = — oo, observing a 
quantity C at t = 0, then proceeding to t = -|-oo, and then going back to 
t = — oo. In the case of the usual T = Green's function with respect to 
the equilibrium ground state, the wavefunction at t = -|-oo is essentially the 
same with the one at t = — oo apart from a phase factor if the initial state 
has no degeneracy [12]. Thus, eq. (PT|) can be decoupled at the time t = +oo. 
However, this simplification does not take place in the nonequilibrium case 
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+ branch 

— oo , 

— oo • 

—branch 

Figure 3: The Keldysh contour for the time evolution. 

of the initial condition eq. ()19p . Therefore, one has to treat the time loop 
including the way back to t — — oo. 

The time-dependent expectation value is defined by using the Heisenberg 
operator Onit), and is written in the form 

{0{t)) ^ Tr[p(0)O^(t)] 

= (f/(-oo,+oo)f/(+oo,t)C»(t)t/(t, -oo))o 

= (f/(-oo,+oo) {T f/(+oo,-oo)a(t)})o 

= (Tc f/,a(r))o. (22) 

Here, the relation among the Schrodinger Os, interaction 0{t), and Heisen- 
berg Onit) representations are given by 

d{t) = e'^^'Ose-'^^' , Onit) = U{0,t) d{t) U{t,0) . (23) 

In eq. fl22|) . Tc and Uc express the time order and time evolution along the 
Keldysh contour, respectively, and t~ denotes the time in the —branch in 
Fig. El 

The perturbation expansion with respect to 7^2 can be carried out by sub- 
stituting eqs. (fT7|) and (fTH|) . respectively, into f/(+oo, — oo) and f/(— oo, +oo) 
in eq. (|22I), as 

m)) = T.T. —T- / <■■■< / dt,...dt^ 

„=0m=0 •'-'^ •'-•^ 

X ( {t n^iA) ■ ■ -7^2(4)} {t n^iti) ■ ■■n2{tm) d{t)} )o . (24) 

The Wick's theorem is applicable to the average (■ ■ ^ because Tii has a 
bilinear form. However, because f/(+oo, — oo) gives a factor (— z)"^ for the 
m-th order terms while f/(— oo, +00) gives a factor for the n-th order 

terms, four types of the Green's functions are necessary to distinguish the 
contributions from these two branches. 




8 



Lif shit z-Pit aevskii 


G— 


G++ 


G+- 


G-+ 


Alternative notation 






G> 


G< 



Table 1: Correspondence between two standard notations. 



2.4 Nonequilibrium Green's function 

We now introduce the four types of the Green's functions, which are required 
in the Feynman-diagrammatic approach to the perturbation expansion along 
the time loop, 



Gtj~^(ti,t2) = 



-*(Tcf/eC,.(tr)4(^2))0 

-i{'fc,^{ti)c]M) 



(25) 
(26) 
(27) 
(28) 



Here Cj^(ti) and c]g.(t2) are Heisenberg operators. denotes the time + 
or —branch in Fig. El For these Green's functions, there is another notation 
used widely in hteratures, and the relation between that and the present one 
by Lifshitz-Pitaevskii [S] is summarized in Table 1. 

Based on the Feynman-diagrammatic approach, the Dyson equation can 
be expressed in a 2 x 2 matrix form; 



Gijiuj) 



rri] 



Im 



G 



G. 
G1 



ij 



■'Im 



''Im 
'Im 



''Im 

^++ 

^Im 



(29) 



(30) 



Here, g^j is the Green's function determined by the unperturbed Hamiltonian 
Til and density matrix p{—oo) for the initial isolated system. The Fourier 
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transform has been carried out for stationary states, 



(31) 



For instance, in the noninteracting case Ti^ 



— 0, the self-energy correction 
is caused only by the couphngs between the sample and reservoirs T^mix, 



,N+lSm,N ) 



Note that the four types Green's functions are not independent, 

G— + G++ = G+- + G~+ , + S++ = - S-+ - 



(32) 



(33) 



Thus, the Dyson equation eq. (j^^ can be expressed in terms of three inde- 
pendent quantities by carrying out a Unitary transformation P~^GP; 





Gl 



1 

V2 



F 



1 1 
-1 1 





9lj 



+ E 

Im 



9ti 

T 

9ii 



ya 



yr 





^ rnj 

F 



(34) 



(35) 



Here, G^ and G"^ are the retarded and advanced Green's functions, respec- 
tively. 



yr 



G- 



G' 



1] 



+ G^ 
+ 



(36) 
(37) 



The function F and Q link closely to a nonequilibrium distribution. Alterna- 
tively, the original four Green's functions can be expressed with these three 
functions, as 



G- 
G- 



[F+(G^ + G'^)]/2, 
[F_(G^-G")]/2, 



G+^ 
G+' 



[F- (G" + G'^)]/2 , (38) 

[F+ (G" -G'')]/2 . (39) 

Similarly, the four self-energies are written in terms of S^', and Q, 

S— = [fi + (S'' + S'')]/2, S++ = [1]- (S" + S")]/2 , (40) 

E-+ = -[^]-(S"-S")]/2, S+- = - (S" - S")]/2 . (41) 



10 



The Dyson equation for the three functions are deduced from eq. (j35|) 

G'- = g'- + / C , = g'' + g" S'' G"" , (42) 

F = F^ + F^T." G'' + g"^ F + g' VlG'' . (43) 

Here, we have suppressed the subscripts for simphcity, and these equations 
should be understood symbohcally. Eq. can be solved formally by using 
eq. (jl2I), as 

= G"{/}~^F°{^"}"^G" + GTIC. (44) 



2.5 Green's function for the initial state 

The unperturbed Green's function Qj^j is determined by Tii and the initial 
density matrix p{—oo). Initially at t — > oo, the three regions are isolated and 
noninteracting. Therefore, g^j.,^ for u = L, R, C is given by 

9ijA^) = 1^ ^ , 9ijA^) = 1^ —T , (45) 

F%A^) = [1 - 2 fA^)\ \9\rA^) - 9trA^) 1 • (46) 



Here, e„;y and (f)n;u{'i) are the one-particle eigenvalue and eigenstate of Hi-^- 
The information about the statistical distribution is contained in the function 
Fij-u via fu{^^) = — fJ^u), where /(e) = [e^^ + 1]"-^. In the system we 
are considering, each of the reservoirs {u = L, R) has a continuous energy 
spectrum, and the isolated sample (z/ = C) has a discrete energy spectrum. 
Thus, the full Green's function becomes to depend only on fiL and fiR, and 
does not depend on fic [El- This is because the contribution of F° to the 
corresponding full one F arises in a sandwiched form {g^}~^ F^ {g""}'^ as 
described in eq. (I44p . Thus, the singular contributions of 6 functions in 
F^ oc [g^ — g"'] of the sample region are canceled out by the zero points of 
the inverse Green's functions in both sides, to yield {g^ }~^ F^ {9"^}^^ = for 
iy = C. 



2.6 Nonequilibrium current for noninteracting electrons 

The nonequilibrium average of the charge and current can be deduced from 
the Green's functions. For instance, by using eqs. and (jSH), an equal- 
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time correlation function can be written in the form 

(4c,J = -^^/^(O.O) = GT+M . (47) 

Therefore, the current flowing from the left lead to the sample is given by 



Jt 



J-oo ZTT 



01 



G 
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(48) 
(49) 



The expectation value for the current flowing from the sample to right lead, 
Jr, can also be written in a similar form. In the noninteracting case Ti^ = 0, 
eq. can be rewritten in terms of retarded and advanced Green's functions 
which link the two different interfaces of the sample ji]; 



IL J — oo 



%{uj) = ATl{uj)GI^{u:)Tr{u:)G 



iVl' 



rL(^) 



Im 



Im 



4^]v+17V+l(^) 



(50) 

(51) 
(52) 



Note that {Jl) = {Jr) (= (J)) in steady states. The outline of the derivation 
are provided in Sec. l3.3l for a single Anderson impurity. Equation implies 
that the current is determined by the electrons with the energy fiR 



<UJ< 



at low temperatures, where fJ^L — IJ'R = ^V- 

For interacting electron systems, the nonequilibrium current can not gen- 
erally be written in the form of eq. (jSUj) . It does only in a particular case 
where the connection of the two leads and the sample has a special sym- 
metry described by a relation r^(e) = Ar^(e) in the notation used in Ref. 
[T^ . In the interacting case, the imaginary part of the self-energy caused 
by the inelastic scattering becomes finite. It links with the contributions 
of the vertex corrections, and the formulation becomes somewhat compli- 
cated. Nevertheless, in the linear-response regime, the dc conductance for 
interacting electrons can be expressed in a Landauer-type form quite gener- 
ally ^1 ^1 even in the case without the special symmetry mentioned 
above [1^]. Specifically, at zero temperature T = 0, the imaginary part of the 
self-energy and vertex corrections for the current become zero at the Fermi 
energy u = 0, and the transmission probability can be written in the form 
of eq. with the interacting Green's functions. We discuss the details of 
these points in Sec. El 
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Figure 4: Anderson impurity connected to two leads 

3 Out-of-equilibrium Anderson model 

We now apply the Keldysh formalism to a single Anderson impurity con- 
nected to two leads as illustrated in Fig. El It corresponds to the = 1 case 
of the Hamiltonian eq. and has been widely used as a model for quantum 
dots. For convenience, we change the label for the sites: the new one for the 
impurity site is given hy Q => and that for the interfaces at the left and 
right leads are ^ L and A^ + 1 ^ -R, respectively. 

3.1 Green's function for the Anderson impurity 

The self-energy in the interacting case can be classified into two parts. 

= S°.(cu) + Sc/(a;)(5,„,5,-rf . (53) 

Here, S° corresponds the one defined in eq. (jH^ . which represents the effects 
purely due to the mixing with reservoirs and sample. The remaining part S;/ 
contains the contributions of the onsite Coulomb interaction U . Substituting 
the self-energy eq. flS^ into the Dyson equation (j^ . we obtain a set of 
equations for the impurity Green's function, 

Gdd{^) = Qddi^) + 9dd{^) Gdd{^) 

-Vl 9ddi^) '^3 Gidiuj) - Vr Qddi^) T3 GRd{uj) , (54) 

Gidi^^) = -VLgL{uj)T^Gdd{i^) , (55) 
GRd{uj) = -VRgR{uj)T^Gdd{^^) ■ (56) 

where ql and qr are the Green's function at the interfaces at left and right, 
respectively. The explicit form of the unperturbed Green's function at the 
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impurity site is given by {Qddi^)} ^ = {^ — ^d)'T3 with one of the Pauh 
matrices T3. Substituting eqs. and into eq. (jSH), we obtain 



ct(uj) 



9dd{^)+9dd{^) ['^M Gdd{^^) 



Therefore, 



{G'Si^)}'' = {9d,i^)r' 



criuj] 



(57) 
(58) 

(59) 
(60) 



Here GfJ{uj) is the Green's function for the noninteracting case. Further- 
more, an alternatively form of the Dyson equation G = g + CE g yields 



Gdhi^^) = -VLGdd{'^)r'igL{uj) 
GdR{uj) = -VRGdd{'^)T'igR{uj) . 



(61) 
(62) 



Thus, the inter-site Green's functions can be deduced from Gddil^) by using 
eqs. and 

The voltage-dependence arises via the unperturbed Green's functions for 
the leads aX v = L and R, 



g:;{oo) 



>-i 



(63) 



F',{u) = [l-2Muj)][g:{uj)-g:{u;) 
= -2t[l-2f,{uj)]r,{uj)/vl, 



(64) 



where T^^u) = — Im [(yf^(ci;)]. The pure mixing part of the self energy cr(u;) 
can be calculated from eqs. and (jH^ . 



^]W(cu) a'-(cu) 
a'^iu) 



(65) 

(66) 
(67) 
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and (t"(u;) = {a'^{u!)}*. Then the noninteracting Green's function G 
can be determined via eq. ()60|) . 



(0), 



dd 



00) 



Gdd 
'-'(id I 



1- /cffM]G, 



(0)a 



where 



(0)a 
dd ' 



(68) 
(69) 
(70) 

(71) 

(72) 
(73) 

and G^^}"'{u}) = ^G^^]\uj)^ . Thus, the effects of the bias voltage arise 
through the distribution function /efr(ci;). 

The full Green's function for the impurity site can be expressed, using 
eqs. (IH-dll and dS!-®, as 

^ (74) 



1 



GU^ 



(75) 

(uj) and 

Fiid{uj) are pure imaginary. The four elements of Gddi'^) can be written, 
using eqs. dSHll-dSni) and as 



where Qd(^) = {^ddl^)}*- Note that VLu{uj) = -S^+(cj) - S^' 



G 



dd 



dd 



'^dd 
^dd 



[UJ 



1-Mu;) G^,,{uj) + U{uj)G':,M 



1 - U{uj) G^,(^) 



-/cfr(^)G'dd(^) - 
-U{uj)[GU^)-GU^)], 
I - U{uj)][G^,M - G'^ddi^)] 



(76) 
(77) 
(78) 
(79) 



where G^^ 
by 



(cij) = —{Gj^ (^)}*) ^-iid /cfr(co') is a correlated distribution defined 

1 



/efr(t^) 



/L(^)rL + /fi(^)r^-^S^'(a;) 
+ - ImS[,(a;) 
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(80) 



With this distribution function, the number of the electrons in the impurity 
site can be written in the form 



(81) 



In the equihbrium case fi = fii = /ir, both /efr(i^) and /efr(i^) coincide with 
the Fermi function f{uj), because of the property eq. 



3.2 Properties of the Green's functions at eV = 

We summarize here the properties of the Keldysh Green's function in the 
limit of the zero-bias voltage eV = 0, at which fiL = fJ^R- In this case, the 
four self-energies can be written in the form 



S[/:cq(^) 



(82) 
(83) 
(84) 
(85) 



Furthermore, in equilibrium, F^deqi^^) and flu^eqi^) are determined by the 
retarded and advanced functions, 



Fdd:eq{^) 
^U:eq{^) 



[l-2/(c.) 
[l-2/(c^) 



C^dd-.cqi^) 
S[/:cq(^) 



" C^dd-.cqi^) 
K:cqi^) 



(86) 
(87) 



Specificality at zero temperature, S^^q(co') and S^~q(ti;) vanish, respec- 
tively, a.t uj > fi and u < fi, because of the Fermi function in eqs. and 
(jHK|l . Similarly, the Green's functions G^^^^{uj) and Gj"^7eq('^) ^Iso vanish at 
u > and u < fi, respectively. Therefore, at the equilibrium ground state, 
the usual T = formalism which yields a single-component Dyson equation 



Gdd:cqi^) 



Gdlcq (^) + GZ.cq (^) Sc7:cq(^) G'dd:eq(^) 



.(0)- 



becomes available. 
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3.3 Current through the Anderson impurity 

The operators for the current flows from the left reservoir to the sample Jl 
and that from the sample to the right reservoir Jr are given by 



Jr = ieY,VR Ic^R^d^ - 4ci?a 



(89) 
(90) 



As mentioned in Sec. 12.61 the expectation values can be expressed in the form 



{Jl) = 2e 
{Jr) = 2e 



°° duj 



vl 



vr 



(91) 
(92) 



The inter-site Green's functions in the right-hand side of eq. ()92j) can be 
expressed in terms of Gddiyj) using eqs. (l ^ - lj^Uj) . and as 



P ' Gr^P 

P-' GdRP 



G%, 
G'm Fm 
GdR 



-vr 

-vr 



GdR FdR 

Thus, using eqs. dSHI), dZS)), and we obtain 

VR [Gd^ - G-^\ = ^ [FdR - Fr^] 



9r G'dd 9RFdd + G% 

Gdd 9r 

Gdd 9r <^dd Fr + Fdd 9r 



(93) 
(94) 



jR 



i.9R - 9r) Fdd - Fp, (G^^ - G;^^) 



G^, [ ATl Ul ~fR)-inu- 2(1 - 2fR) ImS 



Similarly, for the right-hand side of eq. (jHH), we obtain 



(95) 



VL [Glj - G,-+] 



— [FLd — FdL] 



-{gl-9l)Fd, + Fl{G^,,-Gl 



ddJ 



Tl G^, G^, [4Tr ih -fR)+^nu + 2(1 - 2/^) ImS[ 



(96) 
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If the two couplings between the Anderson impurity and leads have a property 
Tl{u!) = XTji^uj) with A being a constant [1^], the expression for current can 
be simplified by taking an average, Fr.x + FbX (j^ . as 

, J. _ rL{JR)+TR{JL) 

= ^ r du [h{u)-fR{u)] [-ImG^,H] . (97) 

Note that (J) = {Jl) = {Jr) for stationary states. 



3.4 



Perturbation expansion with respect to 



So far, we have discussed general properties of the nonequilibrium Green's 
functions. To study how the inter-electron interactions affect the transport 
properties, the self-energy must be calculated with reliable methods. 

Because the noninteracting Green's function, which includes the couplings 
between the leads and the Anderson impurity, have already been obtained 
in eqs. flUH jl - lfTOj) . the remaining task is calculating J]i/{u), for instance, by 



taking G^^J to be the unperturbed Green's function. Then, the interacting 
Green's function Gdd are deduced via the Dyson equation (jSH)). 

The self-energy S^(ct;) can be calculated with the perturbation expan- 
sion with respect to the inter-electron interaction 7i^. For generating the 
perturbation series, it is convenient to introduce an effective action. 



E Jdtdt' r^i{t)Kf,{t,t')riAt') 



-U dt 



-r,\+it)r^nit)v\+it)rii+it) 



(98) 
(99) 



(100) 



Here, /7j.(t) = (^ril_{t), ril^{t) ^ is a two-component Grassmann number cor- 
responding to the — and + branches of the Keldysh contour shown in Fig. 
El The Kernel K^^l is determined by the noninteracting Greens function, 

doj 
2^ 



:^KS(..)e-(-'). 



(101) 
(102) 
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In eq. (|l(J(Jj) . the sign for the interaction along the — branch and that for 
+ branch are different. This correspond to the sign arises in eq. (|24j) . and 
it is determined from which of the time-evolution operators, U{+oo, —oo) 
or U{—oo,+oo), the perturbation terms arise. For S'ex(^?^ v) ^q. 
we introduce an external source of two anticomutating c- numbers j}j{t) 
jl--{t), ji+{t) ) following along the standard procedure pUj . 



(103) 



In this formulation, the Green's functions are generated from a functional 

Z[j], as 



Z[j] = I DriWri e^^^''^'") , 
1 6 S 



m 6Mt) 6j^At') 



(104) 
(105) 



In the noninteracting case, the functional integration can be calculated ana- 
lytically 



Z(°)[0] exp 



:E [dtdt'jiit) G^SMjAt') 



. (106) 



The generating functional Z[j] can be rewritten in the form 

Z[j] = e'<-'^"^) Z(%] . 



(107) 



Here, rj and in the action Sjj{rj\rii) has been replaced by the functional 
derivatives 

5 8 

The perturbation series can be obtained by substituting eq. ()106|) into p07|) 
and then expanding e*'^^ in a power series of Su. 
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3.5 Fermi-liquid behavior at low bias voltages 

The out-of-equilibrium Anderson model has been studied by a number of the- 
oretical approaches. In this section, we discuss briefly the low-bias behavior of 
the Green's function and differential conductance, which have been deduced 
from the Ward identities for the Keldysh formalism |2Ij. In equilibrium and 
linear-response regime, the low-energy properties at max(a;,T) <C can 
be described by the local Fermi liquid theory, where Tk is the Kondo tem- 
perature 1121 . The results deduced from the Ward identities show that the 
nonlinear properties at small bias-voltages eV -C Tk can also be described 
by the local Fermi liquid theory. 

The low-energy behavior of ImS[;(ci;) has been calculated exactly up to 
terms of order uj^, {eVy, and T^, 



ImE[,H = --{Aeq(O)}^ |rTiat(0,0;0, 



X 



[i- L + r) 



(109) 



where Tfj„'.cr'fj{uj,uj';uj',uj) is the vertex function for the causal Green's func- 
tion in the zero-temperature formahsm, and Ae^^u) = —Im G^^^.^^^u) / n . The 
parameter a is deflned by a = (aiFi — aijF/j)/(FL + Tr), where and 
are constants which have been introduced to specify how the bias voltage is 
applied to the equilibrium state. Namely, /il = a^eV^ and fiR = —aReV 
with ai + oiR = 1. 

The real part of the self-energy is generally complicated. However, it is 
simplifled in the electron-hole symmetric case for = —U/2, Tl = Vr, and 
Hl = —fJ'R = eV/2. In this case the spectral wight at the Fermi energy 
becomes Aeq(O) = l/(7rA) with A = F^ + F^, and the low-energy behavior 
of the real part of the self-energy is given by 

ReJ:ij{uj) = (l - z-^) 00 + 0{uj^) , (110) 

(111) 



duj 



uj=0 



Here, the constant Hartree term U/2 is included into the unperturbed part, 
and it set the position of the Kondo peak on the Fermi energy + f//2 = 0. 
Therefore, G^{uj) can be deduced exactly up to terms of order u^, and 
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(el/) 2 from eqs. (ITO and (ITTn|l . 

G'^'M ^ 7 7' (112) 



u + i A + i ■ 



U 



2A(7rA)2 

where the renormahzed parameters are defined by 

A = zA, U = z^T 11,^^ {0,0; 0,0) . (113) 

The order f/^ result ^3] can be reproduced from eq. ()112|) by replacing U 
with the bare Coulomb interaction U and using the perturbation result for 
the renormalization factor z = 1 — {3 — 7i"^/4) + ■ ■ ■, where u = U/{ttA). 

This result shows that in the symmetric case the low-voltage behavior is 
characterized by the two parameters A and U. These two parameters are 
defined with respect to the equilibrium ground state, and the exact Bethe 
ansatz results exist for these parameters. The width of the Kondo resonance 
A decreases with increasing U, and becomes close to A ~ {4:/n)TK for large 
U with the Kondo temperature defined by 

Tk = nA^u/{27r) exp[-n^u/8 + 1/(2m)] . (114) 

The Wilson ratio is usually defined by R = Xs/l, where 7 and Xs are the 
enhancement factors for the T-linear specific heat and spin susceptibility, 
respectively ^T]. Alternatively, it corresponds to the ration of U to A, 

R-l = U/{ttA) . (115) 

The Wilson ratio takes a value R = 1 for U = 0, and it reaches i? = 2 in the 
strong-coupling limit U —>■ 00. 

The nonequilibrium current (J) is calculated by substituting eq. ()112p into 
eq. (jnZI)- Then, the differential conductance dJ/dV are determined exactly 
up to terms of order and (eV)'^, 



dJ_ _ 2^ 
dV ~ ~h 



l + 2{R-iy /vrT\2 _ l + 5(i?-l)^ /eV\^ 
3 \X) 4 VX ^ 



(116) 

Therefore, the nonlinear (eV)'^ term is also scaled by the resonance width A, 
and the coefficient is determined by the parameter i? — 1, or U /{ttA). 
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4 Transport theory based on Kubo formalism 



We have discussed in Sec. 12. 61 that in the noninteracting case the nonequihb- 
rium current can be written in a Landauer-type form, as eq. ()5U|) . The 
similar expression has been derived for interacting electrons in a special 
case, when the couplings between the leads and sample satisfy the condi- 
tion r'^(e) = Ar'^(e) in the notation used in Ref. [T3|. For this condition to 
be held, the interacting sites must be classified into the following two groups: 
one group consists of the sites that are connected directly to both of the two 
leads, and other group consists of the sites that have no direct links (hopping 
matrix elements) to the leads. This condition restricts the application of eq. 
(j5(jp . For instance, if there is an interacting site that is connected to only 
one of the two leads, the condition is not satisfied. Therefore, eq. is 
not applicable to a series of quantum dots as illustrated in Fig. Neverthe- 
less, in the hnear-response regime, the Landauer-type expression of the dc 
conductance, eq. (jl39|) . can be derived quite generally without the condition 
mentioned above jl9j . 

In Sec. 14.11 based on the Kubo formalism, we describe the outline of the 
derivation of eq. ()139|) for interacting electrons. Our proof uses the analytic 
properties of the vertex corrections following along the Eliashberg theory of 
a transport equation for correlated electrons [221 123] ■ The many-body trans- 
mission probability T{e) is given by eq. (jl4U|) . and it is written in terms 
of a three-point correlation function. At zero temperature, the imaginary 
part of the self-energy due to the interaction and the vertex corrections for 
the current become zero at the Fermi energy e = 0. Due to this property, 
the transmission probability at T = is determined by the single-particle 
Green's functions as shown in eq. (|14H) . In Sec. 14.21 the current conservation 
law for the correlation functions is described with the generalized Ward iden- 
tity, which expresses the relation between the self-energy and current vertex. 
In Sec. 14.31 we provide the Lehmann representation of the three-point func- 
tions to carry out the analytic continuation formally. It can be also used 
for nonperturbative calculations of T(e). We apply this formulation to a 
finite Hubbard chain in Sec. 14.41 and show an example of the transmission 
probability T(e) for interacting electrons. 
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4.1 Many-body transmission coefficient T(e) 

We now consider the Hamiltonian Tilot defined in eq. (0) again, which is also 
illustrated in Fig. ^ The dc conductance g can be determined in the Kubo 
formalism, and it corresponds to the a;-linear imaginary part of a current- 
current correlation function Kaa'i^ + ^0"*"); 

g = e hm , (117) 

K^^\ivi)= fdT{T,J^{r)J^,{{)))e'-'\ (118) 

JO 

where a = L or R. The retarded correlation can be calculated via the 
analytic continuation Kaa'i^^ + iO^) = -f^aa'(^^OIij/,-^w+io+' "where z// = 2nl/f3 
is the Matsubara frequency. The current operator Jq, is defined by 

Jl = ^ I] ( cLcoa - cLcia ) , (119) 

Jr = ij^^ni (^N+laCNa - CNaCN+la ) ■ (120) 



Here Jl is the current flowing into the sample from the left lead, and Jr is 
the current flowing out to the right lead from the sample. These currents 
and total charge in the sample satisfy the equation of continuity 

Pc = cl^S<^ ' (121) 



dt 



+ Jr-Jl = 0. (122) 



Owing to this property, the dc conductance g defined in eq. pi7|) does not 
depend on the choice of a and a' P^. Note that K^'aiz) = Kaa'{z) owing 
to the time-reversal symmetry of Ti. 

To calculate the cj-linear imaginary part of Kaa'i^^ + ^0+), we introduce 
the three-point correlation functions of the charge and currents. 



$C;,y(r;ri,r2) = {t^ dpdr) c^M) c]>a{r2)) , (123) 
^L;jAr;n,r,) = {T^JLiT)c^^in)c],^{r,)), (124) 
W(^;n,r2) = {T^JR{T)cj^{n)c],^{T2)), (125) 
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where 5p(j = p^j — {p^)- These three functions can be expressed as functions 
of two Matsubara frequencies iu and ie, 

%-,A^;t,,T2) = ^ E %-,A^e,te + tu)e-'^'^^^-^^e-'^^+''^'-^-^-^ , (126) 

for 7 = C,L, R. We mainly consider the electrons in the central region 
assuming jj' G C. In the right-hand side of eqs. (jl24j) and (jl25j) . there still 
exist the creation and annihilation operators with respect to the leads at 
and + 1 in the operators Jl and Jr. The correlation functions that include 
these two sites as one of the external points can be related to those defined 
with respect to the adjacent sites 1 and A^, by using the properties of the 
Green's function at the two interfaces: 

Go,{z) = -gL{z)v^G^,iz) for l<j<N+l ^^^^^ 
Gj,N+i{z) = -Gj^N{z)vjigji{z) for < j < 

Here gL^z) and gii{z) are the local Green's functions at the interfaces of 
the isolated leads, and + 1, respectively. Using these properties, the 
three-point correlation functions for jj' G C can be expressed as 

^^■jj'{ie, ie + iz/) = ^ Gjj^{ie) Kyj^j^iie, ie + iv) Gj^j'iie + iu) , (128) 
j4jiec 

where ^-y-j^ji includes all the vertex corrections. The corresponding bare 
current vertices are given by 

Ag^.^^,(^£,z£ + zz/) = 6,,,, , (129) 

^L-j4n + ^'^) = ^dis^ + '^^) '^1 j4 ) (130) 

^nldi + = ^iii'>'^y + ^N,ji ^N,n ) (131) 

with 

XL{ie,ie + = -ivUg^iie + iu) - g^iie)] , (132) 

Xniie^ie + iu) = ivl[gji{ie + iu) - g^ie)] , (133) 

We now calculate the tu-linear part of Kaa'{z) taking a and a' to be R 
and L, respectively. Using the three-point correlation functions, the current- 
current correlation function KRiiip) can be expressed as 

Kmiiu) = -^EE'^i(^^'^^ + *^'^) *«;ii(^^'^^ + ^^) • (134) 
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Figure 5: Three analytic regions of $^.11(2;, z \ w). 



Paying attention to the analytic properties of the Green's functions, the 
summation over the Matsubara frequency can be rewritten in a contour- 
integral form. Then, by carrying out the analytic continuation iv —>■ uj + iO^, 
we obtain 



■ /(e + 00)- f{e) ] X^l^ (6, 6 + uj) ^%,{e, e + u) 

M(. . , , A ^[3] . , , A L (135) 



°° de 

-00 2m 
°° de 



UJ 



where /(e) = {e^'^ + 1)^^- The superscript with the bracket, that is, [k] for 
k = 1,2, 3, is introduced specifies the three analytic region of ^R-n{z, z + w) 
and Xl{z, z + w) in the complex 2;-plane. These regions are separated by the 
two lines, Im [z) = and Im {z + w) = 0, as shown in Fig El In each of the 
three regions, ^r-u^z, z + w) corresponds to the analytic function given by 



R]ll 



<^R-u{e + iO+,e + uj + iO+) 
$i?;ii(e-^0+,e + cj + zO+) 
<^'i?;ii(e-«0+,e + cj-iO+) 



(136) 



These analytic properties can be clarified explicitly in the Lehmann repre- 
sentation, eq. (jl58j) . provided in the next subsection. Similarly, the analytic 
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continuation of the bare current vertex Xa{ie,ie + iu) for a = L, R is given 
by 

" AW(e,e + ^) = tSavl[g+{e + uj)-g+{e)] 
X^^K^,^ + ^) = ^s^vl[gt{e + uj)-gM] , (137) 
A[3](e,e + ^) = ts^vl[g~{e + uj)-gM] 
where the factor Sa is defined such that sl = —1 and sr = +1. In this 
section, we distinguish the retarded and advanced Green's functions by the 
label + and — , respectively, in the superscript. In the limit of u; — > 0, the bare 
vertices for k = 1 and 3 vanish as X^^^ (e, e + tu) (x cu. In contrast, for = 2, it 
tends to a finite constant A[^l(e, e) = 2sQ,ra(e) with Ta{e) = — Im [g+(e)]. 
Correspondingly, the asymptotic behavior of ^^^^ (e, e + cu) for small to has 
been investigated by using the Lehmann representation of a four-point vertex 
function [221123, and the result is 



{u for k = 1 

finite for k = 2 , (138) 
LU for k = 3 

for a = L,R. Therefore, taking the — > limit in eq. (jll7|) by using eq. 
()135p for KRhiyJ + ^0"*"), we obtain 

g=^-^ fde Tie) , (139) 



Tie) = 2T,ie)^%,ie,e). (140) 

Thus, the dc conductance is determined by the three-point function for the 
analytic region k = 2. The analytic continuation is performed formally 
by using the Lehmann representation in Sec. 14.31 The result shows that 
$^'^]^(e, e) can be expressed as a Fourier transform, eq. p62j) . of a real-time 
retarded product in eq. ()160p . 

Specifically, at T = the conductance is determined by the value of the 
transmission probability at the Fermi e = 0, and it can be written in the 

form [211231231211, 

r(o) = 4 Tl{o) Gr^(o) r^(o) gUo) • (i4i) 

This is due to the property that the vertex corrections for the current vanishes 
at T = and e = 0, as shown in eq. (jl54j) in Sec. 14.21 Furthermore, the 
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reflection probability is given by 



no) 



2zrL(0)G+(0) = 1 - 2^r^(0)G+^(0) 



(142) 



The current conservation T(0) + 7^(0) = 1 follows from the identity in eq. 
()156|) . Similarly, at zero temperature, the Friedel sum rule for interacting 
electrons is given by 



AiVtot = - log[det5] 

m 

where the S'-matrix is defined by 

l-2zrz.(0)G+(0) 



-2zr,.(0)G+v(0) 
1 - 2zr^(0)G+^(0) 



(143) 



(144) 



-22r^,(0)G+i(0) 

In eq. (I143|) . AA/'tot is the displacement of the total charge defined by 

ieC cr 

+ EE [(4c..) - [clcj,] + EE [(4c..) - {clcj^] , (145) 

where (■ ■ ■)]^ and (■ ■ denote the ground-state average of isolated leads 
determined hj Hl and Hr, respectively. 



4.2 Current Conservation and Ward identity 

The inter-electron interactions generally cause the damping of excitations. 
Therefore, theoretically, the self-energy and vertex corrections must be treated 
consistently with the approaches that conserve the current. In this subsec- 
tion, we discuss the current conservation using a generalized Ward identity. 

The generalized Ward identity can be derived from the equation of con- 
tinuity in the Matsubara form — (d/dr) 5pQ + i Jr — i Jl = jSHj, 

d 

- ^ ^cjjiT; Ti, Ta) + 2 ^Rjj'{r; n, rg) - i ^L-jj'ir; Ti, Ts) 

= Sir - ra) (n, r) - 5(ri - r) (r, r^) . (146) 

It can be expressed by using a. N x N matrix representation for jj' G C with 
the Matsubara frequencies, 

iu ^c{i£, i£ + ii') + i ^R^ie, is + iv) — i ^L{i£, i£ + i^) 

= G{ie) - G{ie + iv) . (147) 
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Here, G{z) = {Gjj'{z)} and ^^{z, z + w) 
version of eq. (jl28|) is given by 



{^^■jji{z, z + w)}. The matrix 



^■y{z, Z + W) 



G{z) \j{z, z + w) G{z + w) 



(148) 



z.z + 



Thus, the identity can also be expressed using Aj{z, z + w) = {A, 
w)}, as 

iu Ac{ie, ie + iu) + i Aji{ie, is + iu) — i A^i^ie, is + iu) 

= {G{te + iu)y^ - {G{ie)y\ (149) 

Furthermore, the Dyson equation for the single-particle Green's function can 
be expressed as 



-til - 



XT/0 

ri-C 



''12 



''21 







^22 A* 








(150) 



(151) 




















vlgR{z) 



(152) 



and S(z) = {Sjj/(z)} is the self-energy due to the inter-electron interactions. 
Therefore, eq. (jl49|) represents a relation between the self-energy and vertex 
functions, and this identity must be satisfied in the conserving approaches. 
Carrying out the analytic continuation of eq. p49|) in the region k = 2, 
ie + iu e + uj + iO~^ and ie e — iO^, and then taking the limit of — 0, 
we obtain 



Agke,e) 



A 



[2], 



-2ImV+.(e) -2ImS+(e) 



(153) 



At T = 0, e = 0, the imaginary part of the self-energy vanishes Im S^(0) = 0, 
and then the current vertices become equal to the bare ones. 



(154) 
(155) 
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Correspondingly, $^|ii(0,0) = G^j^{0)2Tr{0)G%^{0) at T = 0, and then the 
transmission probabihty is given by eq. p41|) . Alternatively, at T = 0, the 
analytic continuation of eq. (jl46j) in region = 2 is written in the form 



G+(0)-G-(0) = G+(0) V+ .(0) - V^i.(0) G-(0) 



(156) 



and the (1, 1) and {N,N) matrix elements represent the optical theorem for 
eqs. (fTiT|l and (ITl^ . 

Particularly, for the single Anderson impurity at = 1, eq. ()147j) becomes 
a single-component equation with respect to the impurity site. In the region 



[21 

/c = 2, it becomes $^ (e, e 



[2], 



G- 



G~^(e) in the limit of 



cij — >■ 0. Furthermore, if the mixing terms have a property Ti{e) = XTji{e), 

[21 

e, e) = — A$p fe, e) follows. Thus, in this case the 



an additional relation $ 



[2], 



dc conductance can be written in the form IT3] , 



9: 



2e' 



single 




-ImG+fel 



(157) 



4.3 Lehmann representation for T(e) 

We now show that the transmission probability T(e) can be expressed in 
terms of a real-time retarded product in eq. ()160p via the Fourier transform 
eq. ()162|) . It shows a direct link between the transmission probability and 
dynamic correlation functions. To prove it, we first of all derive the Lehmann 
representation for ^fi-ii{ie,ie + iu), and then carry out the analytical con- 
tinuation. 

Inserting a complete set of the eigenstates, H\n) = En\n), into eq. p25p 
and using eq. p26|) . we obtain 

^R.ii{ie,ie + = ^ (/|cL|m)(m| jR|n)(n|ci„|/) 

Imn 

^ _{ie + iu + Em- El) {iu + E.m - En) 



~ {ie + En- Ei){ie + iu + E^ - Ei) 

~ {iv + Em- En){ie + En- Ei) _ 
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+^ {^\cia\n){n\jR\m){m\c\jl) 



Imn 



X 



{le + El - En)iiiy + En - Er, 



+ 



'ie + Ei- En){ie + ti^ + Ei - E„ 

p — l3Em 



{le + iu + El - Ejj{iu + En - Ejj y 

where Z = Tre~^^. From eq. (jlSSj) . the analytic continuation to obtain 
$gii(e,e + c^) for k = 1, 2, 3 can be carried out by replacing the imaginary 
frequencies ie and iu by the real ones e and u, respectively, with the in- 
finitesimal imaginary parts shown in eq. (jl36|) . Then the same expressions 
for (e, e + uj) can be derived from the real-time functions 

*S!ii(^;^i,^2) 



e{t-h)e{h-t2) {[{cM , clih)} ,JR{t) 
+ e(t, - 1) d{t - 12) ({ci,(ti) , [cUh) , Mt)] } ) , 



(159) 



e(t - ti) e{ti - h) ({cL(t2) , [ciAh) , Jnit)] } 



0it-t2)e{h-h) {[{c,,{ti) , Cl{t2)} ,JR{t) 

0{t2 - 1) e{t - h) ({cL(t2) , [ciM , JR{t)] } 



(160) 



(161) 

where Jnit) = e*^* J/je~*^*, and 9{t) is the step function. The commutators 
are defined by [A, B] = AB - BA, and {A, B} = AB + BA, as usual. The 
Fourier transform into the real frequencies is given by 

dt dti dt2 e^'^*e*^*i e'"''^^ <l>glii {t; ti, ts) 



27r5(e + cu-e') $gln(e,e 



UJ 



(162) 



For example, a time-ordered function 

F{t-hM) = e{t-ti)e{h-t2){jR{t)c,^{h)cl{t2)) (163) 
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Figure 6: Schematic picture of a finite Hubbard chain. 



is transformed into 

^ -1 e-^^- {l\c\„\m){m\jR\n){n\c^^\l) 

Z {e + uj + E^-Ei + iO+){uj + E^-E^ + iO+) ' 

(164) 

Among the three real-time functions eqs. ()159p - ()161|) . the function for the 
region k = 2, that is, ^^'^^(t; ti, ta) 

in eq. ()160|) determines the transmission 

probabihty T(e) = 2TL{e) $^'^]^(e, e). Because the analytic continuation has 
already been done, the real-time correlation links directly to the transport 
coefficient. This formulation can be used for numerical calculations. 



4.4 Application to a Hubbard chain connected to leads 

In this subsection, we apply the linear-response formulation to a finite Hub- 
bard chain attached to reservoirs, which can be considered as a model for a 
series of quantum dots or atomic wires of nanometer size. A schematic pic- 
ture of the model is shown in Fig. El The Hamiltonian parameters defined in 
eq. ((H) are taken as follows. We take tfj to be the nearest-neighbor hopping 
t, and Uj^j^-j^j^ to be an onsite repulsion U . Specifically, we consider the 
electron- hole symmetric case, at which fi = and ed + U/2 = with being 
the onsite energy. We also assume that the two couDlines are symmetric 
= Fij (= F), a: 




Figure 7: The order terms of (a) self-energy and (b) vertex corrections. 
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Figure 8: Many-body transmission coefficient for = 3 (upper two panels) 
and 4 (lower two panels) is plotted vs e/t for V jt = 0.75 for three values of 
U/{2nt); ( — ) 0.0, (-o-) 0.5, and (-•-) 1.0. The temperature is taken to be 
T/t = 0.0 for (a) and (c), and T/t = 0.2 for (b) and (d). 

To examine the effects of the Coulomb interaction, we calculate the self- 
energy and vertex corrections up to terms of order U"^, the Feynman diagrams 
for which are illustrated in Fig. [7| 19J. These contributions satisfy the gener- 
alized Ward identity eq. (jl53|) that corresponds to the current conservation 
law. In Fig. |H1 the results of T(e) for = 3, 4 are plotted vs e/t for 
T/t = 0.75 for three values of f//(27rt); (— ) 0.0, (-o-) 0.5, and (-•-) 1.0. 
The temperature T/t is taken to be (a) 0.0, (b) 0.2 for = 3 in the upper 
panels, and (c) 0.0, (d) 0.2 for = 4 in the lower panels. 

At low temperatures, there are A^ resonance peaks that have one-to-one 
correspondence to resonant states of the unperturbed system. In addition 
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Figure 9: NRG results of the counductance g and and the number of electons 
rid in triple dots = 3 as functions of at T = 0.0. 

to these resonance peaks, two broad peaks of atomic character appear at 
e ~ ±f//2 for large U. The resonance peaks become sharper with increasing 
U at low temperatures as seen in the panels (a) and (c). However, the height 
of the peaks decreases with increasing U. One exception, which happens for 
odd A^, is the Kondo resonance at the Fermi level e = 0. At this peak the 
transmission probability reaches the unitary-limit value 1.0 for any values of 
U, when the systems have the inversion symmetry Tl = Tji together with 
the electron-hole symmetry j27|En|- The width of the Kondo resonance 
must decrease with increasing A^. For even A^, the transmission probability 
T(e) shows a minimum at e = 0. The characteristic energy scale in this case 
is the width of the valley, which eventually becomes the Mott-Hubbard gap 
in the limit of large A^. The high energy profile of T(e) at |e| ^ 2t in the case 
of A^ = 3 is similar to that for A^ = 4. Namely, the high-energy part shows 
no notable A^ dependence. For U/{27it) = 0.5, the upper and lower Hubbard 
levels at e ~ ±f//2 exist inside the energy region corresponding to the one- 
dimensional band of the width 2t. The two Hubbard levels got outside of this 
energy region for f//2 ^ 2t. At finite temperatures, the resonance peaks at 
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|e| ^ 2t become broad, and the peak height decreases with increasing T. The 
structures of the resonance peaks vanish eventually at higher temperatures, 
and then the even-odd oscillatory behavior disappears fT^. 

Recently, there has been a numerical progress in this subject. The con- 
ductance for = 3 and that of 4 have been calculated away from half-filling 
as functions of with numerical renormalization group (NRG). It has been 
clarified that the conductance shows a typical Kondo behavior as seen in Fig. 
13 Namely, the plateau of the Unitary limit g ~ 2e^//i emerges at the gate 
voltages, ed, corresponding to odd-number occupations of the electrons. In 
contrast, the conductance shows wide minimum when the interacting region 
is occupied by electrons of even numbers |ST| 13^. 



5 Tomonaga-Luttinger Model 

Transport through interacting systems in one dimension has been studied 
extensively for quantum wires, organic conductors, carbon nanotube, etc. In 
this section, we provide a brief introduction to a Tomonaga-Luttinger model 
PBj . to take a quick look at the transport properties of a typical interacting 
system in one dimension. 



5.1 Spin-less fermions in one dimension 

We start with the spin-less fermions described by the Hamiltonian, 



Ho 
Hi 



{Ho)o = - ^) 

k 

^ 51 ^1 '^k+q '^k'-q ^k' ^k 
qkk' 



Ci.C 



k^k 



{CkCk)o 



(165) 
(166) 



In eq. p65|l . the ground-state energy for the noninteracting electrons has 
been subtracted. At low energies, the excitations near the Fermi level play a 
dominant role, so that can be linearlized at the two Fermi points k = ±fc^, 



Ho 



aia 



k^k 



(%«'fc)o 



+ J2M-k-kp)\bibj^-{bib,)o 



(167) 



Here (6^) is the operator for the right-moving (left-moving) particles. 
The summation over k in eq. ()167|1 should be restricted in a range \k\ — 
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Figure 10: Linearized dispersion 



kp < kc with the cut-off momentum kc of the order a band width D ~ 
Vpk^ as illustrated in Fig. I5.ll However, we assume that kc oo, and will 



introduce the cut-off for the momentum transfer q, when it is required 
For low-energy properties, the interactions between the electrons near the 
Fermi level is important. Therefore, the interaction Hamiltonian eq. (jl66|) 
can be simplified by taking only the scattering processes in which all the four 
momentums are close to one of the two Fermi points, that is, /c -|- g ~ ±kp, 
k' -q 



±kp, and k 



±kp, into account; 



Hi 




^k+q ^k' 



-q ^k' ^k 



(168) 



For the scattering process with small momentum transfer g ^ 0, there are 
two types of possibilities for the initial momentums ~ fc' and k ~ —k'. In 
contrast, in the case of the back scattering q ~ ±2kp, the incident momen- 
tums must have the opposite sign k ~ —k'. These scattering processes can 
be described by a simplified Hamiltonian, 



Hi 



7^4 + ^2 + ^1 



(169) 



7^4 



94 V- 



2L 



hg ^k'-q ^k' ^k 



qkk' 



qkk' 



q "k'~ 



q ^k' ^k ! 



(170) 



35 



^2 
Hi 



qkk' qkk' 
q'kk 



(171) 



bk' '^fc + 2L ^'^^+9'+2A:^ ^I'-g'-2fc^ '^fc' ^fc • 

(172) 



In eq. (|172p . q' is a small momentum defined such that q = q' ± 2kp. The 
coupling constants should be taken as g2 — — Vq, and gi ~ V2kp- However, 
in the following we treat these three constants to be independent parameters. 
The momentum-transfer cut-off is introduced for the summation over q and 
q'. The Tomonaga-Luttinger model is defined by 



TYtl = Ho + H4 + 7^2 • 



(173) 



Note that for the spin-less model there should be no distinction between H2 
and Hi 

The interactions H4 and 7^2 can be expressed in terms of the density 
operators pi{p) and p2{p) defined by 



Piip) 



k 



(174) 
(175) 



Here (al.afc)o = d{kp—k) and {h\h^)Q = 9{kp+k), and these terms are required 
to define the deviation from the noninteracting value without an ambiguity 
caused by the occupation of the negative energy states. Equations ()17()j) and 
()17H) can be rewritten in the following forms apart from a renormalization 
of the chemical potential that can be absorbed in to kp, 



^4 = y 51 Pii-p)PiiP) + y II P2i-P)p2ip) , 

p>0 ^ p>0 

92 



H2 



r YPli-P)p2iP) + J- Y P2i-P) PliP) ■ 
^ p>0 ^ p>0 



These two density operators satisfy the commutation relations 

Lp 



Piip), Pii-P') 



Lp_. 

Oppi , 
Ztx 



P2{-P) . P2iP') 



271 



(176) 
(177) 



;i78) 
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One notable feature is that these two commutation relations are equivalent 
to those of the bose operators, 



2tt 



Spp' , 



cl = pA-p) 

dp = P2{P) , 



Lp 



Spp' , 



(179) 
(180) 
(181) 



where p > 0. The commutation relation of the density operators and Hamil- 
tonian can be calculated by using eqs. (jl67|) and (jl74|) - (jl77|) as 



Piip) , ^0 = Vppp^{p) 



P2 (-P) , ^0 



VfPP2{-p) , (182) 



Pi(p), = QiVppp^ip) , [p2{-p), H4J = QAVpPPii-p) ' (183) 
Pi(p), ^2] = 92VfPP2{p) > [p2{-p)^ ^2] = 92VfPPi{-p) , (184) 
where = g^/ {2nvp) , and ^2 = 92/ {^nvp) . 

5.2 Two conservation laws 

The operator for the charge and current are defined by 

Pcip) = PiijP) + P2{P) , pAp) = pAp) - P2{P) ■ (185) 

In the real space, the operators for the left and right movers z/ = 1, 2 are 
written in the form 



PvK.x) 



^E(p.(p)e^"^ + p.(-p)e-^"^) 



(186) 



The equation of motion for Pc(^) ^^"l Pji^) derived from the Heisenberg 
equation using the commutation relations in eqs. p82j) - p84|l . 

d d 

— p^{x,t)+vj — pj{x,t) = 0, vj = vp{l + g4- g2) , (187) 
d d 

— pj{x,t) + VN-g^Pc{x,t) = 0, Vn = vp{l + gi + g2) . (188) 
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Because there are two independent equations for p^{x,t) and pj{x,t), the 
exphcit form of these Heisenberg operators can be calculated analytically; 



92 



p^{x,t) = 0, — 



dx^ 



pj{x,t) = 0, (189) 



(190) 



The relation among the three velocities f j, and v can be summarized as 



KpVp , 



^ K 



1+94 + 92 



(191) 



5.3 Charge and current correlation functions 

Owing to the property shown in eq. ()189p . the correlation functions for the 
density operators 

xUiP^t) = ileit){[p^{p,t),p^i~p)]) , for p,u=l,2 (192) 

can also be calculated exactly. The equation of motion for these correlations 
are given by 



L 



^(t) ( PuiP) , Pui-P) 



Bit) 



L 



dpf,{p,t) 



dt 



Pui-P) 



^rljit) + ie{t)^{[[p^{p,t),nT^^] ,pA-p)]) 



2ti 

where is a Pauli matrix: 



(193) 



(194) 



The commutation relation [p^{p,t) , T^tl] in eq. ()193p can be calculated by 
using eqs. p82|) - ()184|) . Then, by carrying out the Fourier transform with 



"0 


r 




'0 -i' 




"1 


0' 




'1 


0" 


.1 


0. 




J 0. 




.0 


-1. 


, 1 = 


.0 


1. 
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respect to t, we obtain 



iU T — V p 



^) = - 1^ 1 , (195) 



COSh(y9 = 



1+^4 



1 / 1 



+ K, 



p ' 



(196) 



sinh if 



92 



1 / 1 



2 _n2 2\Ka 



(197) 



Note that, cosh ipl + sinh = exp (</? r^). The Bogohubov transformation 
given by exp {ip t^/2) has a property, 



3 (lPT^\ n f^T^ 

T = exp T exp 



X 2 ; V 2 

Therefore, eq. ()195|1 can be diagonahzed, as 



(198) 



exp ( \ {uJT^ - Vppl} exp 



With this transformation by exp (y? t^/2), the operators Cp and Clp are 
transformed into 

'cosh(y9/2) sinh((y9/2) 

(200) 





Cp 







cosh((/3/2) 



sinh(y9/2) 



where the bose statistics is preserved for the new operators. 
The exphcit form of x'"(p, cj) is determined by eq. ()199p . 



1p ) 1p' 



(201) 



6pp' . 



X'{p,u) = -Aexpf-^] iDlip,u)T' + DL{p,u;)l\expL^^ 
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= - ^ I L'+(p,cj) + DLip^uj) (coshv^l - sinhy^r^) | , 

(202) 

Dl{p,uj) = i f i ± i -] . (203) 

The charge susceptibihty x^dVi ^)-> which corresponds to the p^-p^ correlation 
function, is given by 

= T.x;.{p,uj) = -— r^4?^V^- ^204) 

Then, the uniform charge susceptibihty is given by limxlip, 0) = Kp/inv ). 

It becomes twice as large for the spin 1/2 fermions. 

The current operator is determined by eqs. p85p - ()187|) as 

J = evjpj. (205) 

Therefore, the J-J correlation function is given by 

K'i,,.) ^.tl^ (206) 

Then, the conductivity can be calculated with the Kubo formula, 

K^{p,uj)-K^{p, 0) 



a{p,uj) 



lUJ 



TT (u; + iSy - {v p) 



2 



(207) 



The conductivity a{p, u) for a uniform p = and stationary u = field 
depends on the order of taking the limits of p — >^ and uj ^ 0. The Drude 
weight corresponds to the p — > limit, 

Rea(0,w) = e^KpV^5{uj) , (208) 

In the real space, the conductivity takes the form 

a{x,u;) = r^a{p,u)e^v^ 

J-oo ZTT 
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dp i e^KpV 



27r 2tt 



+ 



ijj — VpP + z5 + f + ib 



2ti 



(209) 



The dc conductance corresponds to the u; — hmit, 



a(x,0) : 

where has been reinserted. 



27,%^' 



(210) 



5.4 Boson representation of the Hamiltonian 



We have seen in the above that the bosonic excitations play an important 
role on the transport properties of the Tomonaga-Luttinger model. Corre- 
spondingly, there is one notable feature in the commutation relations for the 
density operators in eqs. ()182j) and p83|) : the two parts of the Hamiltonian 
Tio and Tii/gi show the same commutation relations. Therefore, one can 
introduce an effective Hamiltonian Tio defined by 



2-itvt 



p>0 



-p)pl{p) + 



27rt>, 



EP2( 

p>0 



'P) P2{P) 



(211) 



which reproduces the commutation relation eq. ()182j) . Thus, the correlation 
functions can be calculated exactly by using TIq as a replacement for TIq. The 
effective Hamiltonian is written in a bilinear form with the boson operators. 



7^0 + "^4 + 

E [C'p C'_p VpP cosh ifl + sinh ip 

p>0 

E ( 7p + 7-p 7-p ) + const . 



a 



C. 



p>0 



(212) 



In this section, we have discussed only the two-particle correlation func- 
tions. The equation of motion for the single-particle Green's function can 
also be written in a closed form |M1 , and the precise calculations have 
been reported in Refs. [HUl 13 7j . Furthermore, for comprehensive description 
of the bosonization, see Ref. 
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